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Flow  between  concentric  rotating  cylinders  - 
-  a  note  on  the  narrow  gap  approximation 

by 

F.  E.  Bisshopp 

Derivations  of  conditions  under  which  steady  flow  between  con¬ 
centric  rotating  cylinders  becomes  unstable  to  perturbations 
of  infinitesimal  amplitude  are  often  carried  out  within  the 
framework  of  the  narrow  gap  approximation  in  which  the  separa¬ 
tion  of  the  two  cylinders  is  presumed  to  be  much  smaller  than 
either  radius.  Now  a  narrow  gap  limit  in  which  the  ratio  of 
the  separation  to  a  characteristic  radius  approaches  zero  is  not 
uniquely  defined,  but  rather  the  form  of  the  resulting  govern¬ 
ing  equations  depends  upon  the  manner  in  which  the  passage 
to -the  limit  is  made.  A  systematic  survey  of  narrow  gap  limits 
may  therefore  be  useful. 

Let  there  be  a  viscous  incompressible  fluid  (density, 
p  ,  kinematic  viscosity, v)  confined  by  concentric  cylinders  of 
infinite  length  (radii,  R^  and  Rg).  Motion  of  the  fluid  is  to 
be  caused  by  rotating  the  cylinders  at  angular  velocities  Q^and 
Qg,  by  translating  them  in  the  direction  of  their  common  axis 
at  velocities  W-^  and  Wg,  and  by  impressing  external  pressure 
gradients  Pg  =  (3p*/r*0<p*)e  and  =  (8p*/0s+)e  (where  the 
asterisk  denotes  the  usual,  dimensional  variable  and  r*,  cp*, 
z*  are  cylindrical  coordinates  of  a  fixed  system  whose  z*-axis 
coincides  with  the  axis  of  the  confining  cylinders).  We  shall 
choose  a  characteristic  radius  R  which  1.3  of  the  order  R-^  and 
Rg  (we  assume  R^Rg  /  0)  and  define  a  gap  width  parameter  as 

The  results  presented  in  this report  were  obtained  in  the 
course  of  research  sponsored  by  the  Office  of  Naval  Research 
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e  s  g  where  d  =  R2“R1 

A  narrow  gap  limit  Is  then  a  limit  where  e — >0. 


Now  let  us  Introduce  a  set  of  dimensionless  variables 


defined  by 


n2  n-  r2  ni+ 

r*  =  R(l+ex),  cp*  =  e  y,  z*  =  Re  °z,  t*  =  — e  t  , 

„  m,  41  M0  nu  M,  m, 

u*  =  ^  e  Xu,  v*  =  ■gCe  2V+e  2v>,  w*  =  ^(e  ^W-e  ’w)  ,  |(2) 

2  IV  m^ 

p*  =  (e  +e^p)  ,  (n„  >0  and  n.  >  0)  J 

R^  2  J 

vAiere  V(x),  W(x)  and  P(x,y,z)  describe  the  equilibrium  flow. 

The  equations  governing  the  equilibrium  flow  are  obtained  by 
substituting  the  relevant  parts  of  equations  (2)  in  the  Navier- 
Stokes  equations  and  the  equation  of  continuity  expressed  in 
cylindrical  coordinates.  Certain  terms  in  the  resulting  equa¬ 
tions  are  dominated  by  others  when  e — >0  regardless  of  how 
and  n^  are  chosen;  when  such  terms  are^delted,  the  governing 
equations  ares 


v2  . 

II  H^-Mg-ng 


v"  -  £  V2'Ms'n2py  -  s-5 ,2~%; 

-*  .  e v2-«3-3p  .  aL  e2-v; 

z  pv2  J 

where  subscripts  x,y,z  denote  partial  differentiation. 
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An  exponent  such  as  ^  which  can  be  eliminated  from 
all  but  one  equation  can  immediately  be  eliminated  altogether; 
in  this  case  the  requirement  that  V  and  P  shall  be  present  in 
a  narrow  gap  limit  implies  that  the  behavior  of  P  in  the  limit 
where  e  — >  0  is  determined  by  the  behavior  of  V  according  to 
the  relations, 


\  *  2M2  +  1  , 


Px  =  V 

Solutions  of  equations  (3)  for  V  and  W  under  the  impressed 
boundary  conditions  are 

V  -  6  *£“  (Q+xAQ)+  |  ^5  PgCx-x^Cx-xO  j 

P  v 


~M  2 

W  =  e  3C~(WfxAW)+  |  ^2  P5(x-x1)(x-x2)j  , 

pv 

where  x^  and  Xg  are  defined  by 


r1„(x1)  =  R^  and  r^Cxg)  =  R2  > 

i.e. 

x^  =  -(R-R^)/d  and  Xg  =  (R2~R)/d  , 

and  where 


Q  =  x2Ql-x^Q2,  AQ  s 

(7) 

W  =  x2W1-x1W2 ,  AW  =  W2-Wx  . 


The  nonlinear  equations  governing  time -dependent 
perturbations  of  the  steady  flow  again  include  terms  which  al¬ 
ways  vanish  in  the  limit  where  e — >  0;  after  such  terms  are 
deleted  the  equations  become 

m-,-1  m?-n?  nn-n., 

e  u v;  +  e  v  +e  J  J  w_  =  0  , 
a  y  Zt 


m. 


“V 


m. 


M0  2  m0-n„ 


e  Lu  =  -  e  px+  (e  v+  e  V)  -2e 

m0  Mo+m-,-1  t  m^-n5  m-,-n9 

e  ^Lv+e  *  V  u  =  -  e  Pv+2e  uv 

v  i7 


and 


nu  M,+nu  -1  t  m^ -n., 

;  JLw+e  J  W  u  =  -  e  Ji 


(8) 


where  the  prime  denotes  , 


a.2  «  ♦  E"2'n2vfc  ♦  e”3"’34;  , 


and 


2  _  -2  ET_  J 
V  "  e  2  +  £ 

ax^ 


~2n2  a2 
ay2 


ay  1  e 

-2n,  ? 


3z 


y  C9) 


+  e 


az 


The  linearized  equations  governing  an  infinitesimal  perturbation 
of  the  form 


u  a  u(x)exp  (i(k2y+k^z)+\t  ]  , 
p  *  p(x)exp  Ci(k2y+k^z)+\t  j  , 


(10) 


are 
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m,  —1 


JU-I  —  j.  LJnlirt  U1-  — ll* 

e  Du+i(e  d  ^k^H-  e  *  ^k^w)  =  0  , 


Lu 


m„-ru  m-j-n 

c 2v+  e  *  ^ 

ml+-m1“l 


e 


Dp  +  2  e 


M0+m0-m,  m0-m,  -n0 

d  *  vv-21  e  k^v, 


nv -mo-n0  M^+m-,  -m0-l  t  m,  -nu-n5 

Lv  =-ie  *  ^k2p-  e  V  u+2i  e  A  *  ^kgU, 


and 


m, -m-j-n,  M^tm-i  — m.,“l  ( 

Lw  =-i  e  ^k^p  **'  e  ^  ^  W  u, 


(ID 


where  D  =  ^ 


and 


L  = 


— nj,  Mp-n?  M^-n^ 

5  \+  i(e  k0V+  e  3  ^-.W)  + 


+  u'^l*  e-V) 


■J 


(12) 


In  addition  to  equations  (11)  the  velocity  components  must 
satisfy  the  boundary  conditions 


u  =  v  =  w  =  0  at  x=Xj  and  x=x2  (13) 

With  regard  to  choice  of  the  exponents  m^-l, 
m2_n2 ’  8111(1  m3“n3>  we  see  immediately  that  at  least  two  of  them 
must  be  equal,  for,  otherwise,  the  equation  of  continuity  implies 
that  the  only  solution  is  the  trivial  one  where  u=v=w=p=0. 

The  case  where  two  of  the  above  exponents  are  equal  and  greater 
than  the  third  simply  cannot  occur,  for  then  the  equation  of 
continuity  implies  that  the  veolcity  component  associated  with 
the  third  exponent  is  identically  zero;  this  means,  in  effect, 
that  we  have  tried  to  chose  the  exponent  bl  ;  associated  with  the 
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third  component  too  small.  Now  the  exponents  m^  indicate  the 
behavior  of  the  leading  terms  in  expansions  of  the  pressure 
and  the  velocity  components  in  powers  of  e ,  and  consequently  a 
value  for  which  is  too  small  in  the  above  sense  is  not  allow¬ 
ed  since  the  leading  term  of  an  expansion  is  by  universal  agree¬ 
ment  not  identically  zero.  Thus  a  rough  classification  of 
possible  limits  on  the  basis  of  the  continuity  equation  contains 
four  cases*  the  case  n^-l^mg-ng^m^-nj  and  the  three  cases  where 
two  of  the  exponents  are  equal  and  less  than  the  third.  We 
shall  first  treat  the  case  where  the  three  exponents  are  equal; 
then  it  will  be  shown  that  the  remaining  cases  do  not  occur. 

If  we  now  take  m^-lamg-rig^m^-n^  a  (say),  equations 

(11)  become 


Du  +  KkgV+k^w)  =  0, 


and 


mr-a-2  il,+n9-l 

Lu  =-e  Dp  +  2  e  *  *  Vv, 

A 

ift  -a-2n0  M~-n0  , 

Lv  =»i  e  v  TtgP  -  e  *  v  u  , 


mu -a— 2n-,  M^— n^  , 

Lw  =  -i  e  ^  ^k^p  -  e  J  %  u  , 


(14) 


where  two  terms  have  been  dropped  since  mi“m2-n2!=~2n2+‘l‘  and 
m2 -mi -n 2= -1.  As  we  shall  see  presently,  the  exponent  a,  which 
now  remains  unspecified,  can  be  assigned  if  one  considers 
secondary  motions  of  finite  amplitude. 
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The  exponents  we  have  to  deal  with  are  Mg-n2, 

Mg+ng-lj  -2n2,  -2^,  -n^,  and  n^-a.  Of  these,  the  last 
two  are  determined  by  the  others  In  virtue  of  the  following 
considerations!  (lj-n^  appears  but  once  (in  L)  and  therefore 
must  be  set  equal  to  the  smallest  of  the  remaining  exponents  in 
equations  (14)  in  order  that  the  parameter  X  shall  appear  in 

the  limit  and,  (2)  m^-a  appears  only  in  conjunction  with  p  and 

therefore,  in  order  to  insure  that  p  is  the  leading  term  of  an 

expansion  in  powers  of  e,  must  be  given  the  minimum  value  con¬ 

sistent  with  the  requirement  that  p  shall  not  vanish  identically. 
Furthermore,  we  must  add  the  restriction  that  none  of  the  ex¬ 
ponents  which  appear  in  L  shall  be  less  than  -2;  otherwise, 
equations  (14)  can  be  reduced  to  a  single  second-order  dif¬ 
ferential  equation  for  u  which  can  be  shown  to  have  no  nontrivial 
solution  which  satisfies  the  boundary  conditions  u=Du=0  it 
and  Xg.  Thus  we  have 

n2  <1  ,  n~£l  ,  Mg^.ng-2,  M^>.n^-2,  nlf=2.  (1?) 

If  we  now  add  to  the  above  restrictions  the  requirement  that 
the  functions  V  and  W  shall  in  fact  appear  in  the  limiting 
equations,  we  obtain  the  relations, 


and 


or 


M^-n^  =  -2  with  n^  <  1 


Mg-n2  =  -2  with  2  —  n2  - 1  y 


Mg+ng  ■  -1  with  n2  » 


(16) 


8 


from  which  it  follows  that  a  choice  of  n2  and  n^  completely 
determines  the  remaining  exponents. 

Now  let  us  consider  the  case  where,  x^n^l*  If  we 
make  the  transformation 


w  — > 

w/k^  , 

v  — >  v/k2  , 

W  — > 

W/k3  , 

V  -»  V/k2  , 

(17) 


(20) 
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with  the  boundary  conditions 
u  =  Du  =  0  at 


x  =  Xj  and  x  =  x2 


(21) 


If  we  now  consider  the  case  where  n^  <  1  we  obtain 

governing  equations  which  have  the  same  form  as  equations  (18) 

through  (21)  except  for  a  replacement  of  the  quantity  k7  by 
i-n  3 

e  ->  throughout.  Thus  it  is  clear  that  no  new  narrow  gap 

limits  can  be  obtained  by  taking  n^  <1  since  they  are  all 
covered  under  the  case  n^=lj  k^=0.  By  means  of  a  similar 
argument  we  can  establish  that  the  cases  where  ^  <n2  <1  are 
equivalent  to  the  case  where  n2»l,  k2=0* 

When  n2=^  a  new  term  is  introduced  in  the  govern¬ 
ing  equation  for  u,  and  a  second  narrow  gap  limit  is  obtained* 
The  governing  equations  then  become 


Du  +  i(v+wj  =  0 
Lu  -  -Dp  +  2VvA| 


Lv 


-V  u 


Lw  =  -ik2p  -  W*u 


(22) 


where 


3 2  2  ^ 

V  *  e2k2  {“  (&-xaq)+  §  P2(x-x1)(x-x2)j 

W  =  ek3  {^(WtxAW)+  \  £^P3(x-x1)(x-x2)3 


pv 


L=\+  i(V+W)+  (k2-D2) 


(23) 
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The  pressure  and  a  velocity  component  may  be  eliminated  from 
equations  (22)  to  give  the  equations 


(LCkf 


■•D2) 


,  it  n 
i(V  +W 


k-  2 

)\  u  =  2(r^)  Vu 

^  A>a 


Lv  =  -  v’u 

with  the  boundary  conditions 


(24) 


u  =  Du  =  v  =  0  at  x=x^  and  x-Xg  (’25’) 

When  n2  <  ■£,  the  exponent  M2-n2  must  be  greater 
than  -2  and  we  then  obtain  the  equation 

(X+iW+(k2-D2) )v  =  0  (26) 

From  the  equation  above,  its  complex  conjugate  and  the  boundary 
conditions  v(x1)=v(x2):a0 ,  we  obtain  the  relation, 

X+X*  -  -f  (k2  |v  |2+ |Dv  |2)dx/J  |v|2dx,  (27) 

°X1  /  X1 

from  which  it  follows  that  there  can  be  no  unstable  solution 
unless  v®-0.  Thus  we  find  only  the  two  limits  given  above  if  we 
require  that  all  three  velocity  components  shall  have  nontrivial 
expansions  in  powers  of  e  with  leading  terms  which  appear  in  the 
limiting  form  of  the  continuity  equation. 

The  exponent  a  which  has  been  left  unspecified 
may  be  assigned  by  requiring  that  equations  (8)  shall  include 
nonlinear  terms  which  limit  the  exponential  growth  of  unstable 
perturbations  predicted  by  the  stability  equations.  It  is 
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readily  verified  that  a  *=-2  Is  the  only  choice  possible. 

Narrow  gap  limits  in  which  one  of  the  velocity 
components  does  not  appear  in  the  lowest  order  continuity  equa¬ 
tion  fall  into  two  categories t  (1)  limits  which  can  be  obtained 
from  the  two  limits  we  have  discussed  above  by  setting  one  of 
the  velocity  components  to  zero  and  (2)  limits  which  introduce 
new  terms  in  the  governing  equations.  It  can  be  shown  that 
there  are  no  nontrivial  solutions  of  the  equations  governing 
the  limits  in  the  first  category,  and  thus  we  need  consider 
only  the  seoond. 

From  equations  (11)  it  follows  that  no  new  terms  can 
be  introduced  into  the  governing  equations  when 


(m^-n^)  “  p  >(m2-n2)  ■  (mj-nj)  »  -2; 

hence  we  need  not  consider  this  case* 

When 

(n^-l)  «p>(m2-n2)  -  (m^-n^)  *  "2  , 

we  have 

Mg+m^-Og-l  B  M2^i2+(p*fr2)  ^  M2*^2 

M^+m^-fflj-l  «  M^Hi^+(p+2)>  M^-n^  r 
m^-m2-n2  *  (p+2)+l-2n2  >  -2n2 


(28) 


(29) 


(30) 


In  Virtcp  of  equations  (29)  and  (30)  the  governing  equations  are 
reduced  to 
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If  the  above  system  is  not  to  degenerate  to  a  lower  order 
system  (for  which  the  boundary  conditions  cannot  all  be  satis¬ 
fied)  ,  the  remaining  exponents  must  be  chosen  in  such  manner 
that  no  power  of  e  appearing  in  the  above  equations  is 
less  than  -2.  Furthermore,  if  V  is  to  affect  the  limit  at  all 
we  have 


Mg+mg+l  =  -2  ,  (38) 

and  if  the  third  equation  is  not  to  become  Lv  =  0  we  have 
either 


n^-m2-n2  =  -2  , 

(a) 

Ml“m2”^  “  —2  , 

Cb) 

-mg-hg-l  ■  -2  « 

(c) 

(39) 
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It  can  be  shown  that  equations  (38)  and  (39)  (a,b,  or  c) 
cannot  be  satisfied  when  the  indices  satisfy  the  inequalities 
implied  in  equations  (37)  and  (34).  Thus  there  is  no  narrow 
gap  limit  here,  and  this  exhausts  the  possible  cases. 

The  two  narrow  gap  limits  are* 

Limit  It  n2*n^=l,  n^=2,  n>4=-2, 


m2=ro3=M2=M3=-l  and  M^=-l. 
The  governing  equations  are 


and 


where 


jL(k^+k3-D2)+i(k2Vn+k3w")ju  =  0 


u  «  Du  =  0  at  x-x^  and  x?=x2  r 
L  *  X  +i(k2V+k3W)  +  (k2+k2-D2) 


(40) 


There  are  the  following  six  Reynolds  numbers  to  be  considered 
in  general* 


RSd  RAQd  1  d3  r  Wd  AWd  nnj1  Id3  r 

v*  v  *  2  2  r2*  *  oiia  2  2*3 

pv  v  v  pv 

The  first  and  fourth  of  the  above  may  be  eliminated  by  a 


(41) 


redefinition  of  X. 

Limit  2)  n2=2»  n3=1»  nlf“2»  nVs“2» 


m3“M3”“1  ^““2 
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The  governing  equations  are 

[L(k2-D2)  +  l(k2v"+k3wn)]u  »  2k2Vv, 


and 


where 


Lv  =  -  V  u. 


u=v=Du=0  at  x=x-^  and  x^Xg  y 


L  =  X+i(k2V+k3W)+(k2-D2) 


(42) 


Here  the  six  Reynolds  numbers  are* 

(43) 

In  view  of  the  fact  that  the  Reynolds  numbers  as¬ 
sociated  with  V  approach  zero  for  fixed  Q ,  AQ  and  P2  faster 
in  limit  2  than  in  limit  1  by  the  factor  / d/R ,  it  might  be 
surmised  that  limit  1  is  the  more  important  of  the  two.  In 
fact  this  would  be  exactly  the  wrong  conclusion.  Let  us  con¬ 
sider,  in  order,  the  three  cases  where  W  vanishes  and  the 
dominant  term  in  V  is  due  either  to  Q,  AQ  or  Pg.  In  the 
first  case  the  fact  that  the  Reynolds  number  associated  with 
Q  can  be  eliminated  in  limit  1,  but  not  in  limit  2,  indicates 
that  it  is  only  limit  2  which  has  any  bearing  on  the  problem. 

lO  p 

In  the  second  case,  the  neglect  of  2  Rd'Pg/pv  relative  to 

p 

R  AQ/v  in  limit  1  gives  the  problem  of  pure  shear  flow  between 
parallel  planes.  Since  this  problem  is  known  to  have  only 


v 


d 

R 


V 


d 

\J  R* 


1 

2 


pv 


2  *2  \ 


d 

R’ 


Wd  AWd 


and  \ 

cl  (. 

pv 
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stable  solutions,  it  is  again  limit  2  which  matters.  It  is 

only  in  the  third  case  that  limit  1  actually  prevails; 

♦ 

however,  as  Lin.  has  already  remarked  ,  if  one  adopts  the 

narrow  gap  approximation,  wherein  the  equations  of  the  narrow 

gap  limit  are  used  for  finite  but  small  e ,  then  on  the  basis 

of  the  computed  critical  Reynolds  numbers  for  the  two  limits 

it  is  found  that  limit  1  is  dominant  only  when  e  is  less  than 
| 

0.26  x  10  ;  and  this  is  indeed  a  very  narrow  gap. 


* 

C. 


C .  Lin  Theory  of 


Stabllll 


